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LOWER BOUNDS FOR THE BUCKLING 

* 

PRESSURE OF SPHERICAL SHELIA 
By Nicholas J. Hoff and Tsai-Chen Soong 
Stanford University 


SUMMARY 

/ <? 

An approximate large-deflection theory is developed for the calcula- 
tion of the hydrostatic pressure under which a thin-walled spherical 
‘ shell buckles. The theory is based on the assumptions that the buckle 
has the shape of a spherical cap and that a yield hinge develops along 
the circumference of the cap. The effect of the test arrangement on the 
buckling pressure is taken into account. In spite of the rough approx 
mate nature of the analysis reasonable numerical values are obtained for 
the buckling pressure. 

V 

INTRODUCTION 


In a review of the state of knowledge regarding the buckling of 
thin elastic shells, Fung and Sechler^ ** noted that a complete 
spherical shell buckles in the form of a small inward dimple of small 
solid angle at some point on the surface of the shell; a fact entirely 
at variance with the linear theory of buckling". Bie non-linear theories 
developed to explain the discrepancies have all dealt with shallow 
spherical caps rather than with complete spheres and almost all the ex- 
perimental effort related to spherical shells was made with spherical 
caps supported along their edge. The only known test with a complete 
spherical shell was carried out by Sechler and Bollay and was reported 
by Tsien. ^ In this test a brass shell buckled under approximately one- 
quarter the critical value of the external pressure of the linear theory. 
More recently, Thompson^ published the results of two experiments carried 
out with polyvinyl chloride shells of a radius to thickness ratio of 
approximately 20. Buckling occurred after about three-quarters of the 
classical value of the buckling pressure was reached. The paper also 
contains a large -deflect ion strain energy analysis which yielded satis- 
factory agreement with the results of the experiments. 

*The work here described was carried out under NASA Grant NsG-93-60. 

^Superscript numb ers indicate entries under References at the end of 
the paper. 
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In a recent paper 4 the senior author developed a rough approximate 
theory of the buckling of axially compressed thin-walled c ^ c ^® r 1 ' 
drical shells on the basis of the assumption that plastic hinges f 
S the ridges when the shell snaps into the large-deflected shape. 
Since reasonable agreement was found betweenavailable e^erimentalda 
and the results of the approximate theory, it appea 

apply the same method of analysis to the spherical shell. The a ° 8an ^ 
o? noticeable permanent deformations after the pressure is released ^ e _ B 
not prove that yielding had not taken place in 4he . ■, the 

ficient amount of elastic energy is stored in the buckled shell, 
shell can snap back into a shape that appears to be the same as the 
initial one. Nevertheless some damage is always done because the buck- 
ling experiment cannot be repeated many times with the same specimen. 

Since the details of the loading machine were found to 
significantly the buckling load in the case of the cylindrical shell, 
the loading apparatus to be used with the spherical shells is carefully 
deNnedl^the present paper. The numerical values obtained for a 
special case and plotted in two figures appear to be reasonable. Never- 
theless no claim is made for the accuracy of the theory because it is 
based on rather far-fetched assumptions. It obviously cannot yiel 
than approximate lower bounds on the buckling pressure because it con- 
sists of equating energy quantities before and after buckling. Th 
relative success of the present approach may, however, indicate that the 
yield stress is one of the physical factors to be taken into account in 
more rigorous analyses of the large-deflection buckling process. It 
was included in the empirical formula proposed by Donnell tventy-eigh 
years ago. 


SYMBOLS 


C correction factor for influence coefficient (see Eq. 23) 

d radial components of displacement 

h wall thickness of specimen 

H wall thickness of container 

E Young' s modulus 

k influence coefficient (see Eqs. 21 and 22) 

K Baker- Cline multiplier (see Eq. 23) 
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moment resultant 


edge moment resultant 
pressure 

radial component of force resultant 

non-dimensional quantity defined in Eq. 24 

radius of specimen 

radius of container 

volume inside specimen 

volume between spheres 

work 

non-dimensional quantity defined in Eq* 25 

semi-vertex angle of bulge 

multiplier defined in Eq. 62 

ratio of specific heats of air (= 1.4) 

lack of fit defined in Eq. 20 

curvature 

Poisson 1 s ratio 

stress 

rotation 

In addition, the following subscripts are used: 

atmospheric 

adiabatic 

inside specimen after loading; also refers to cap 
inside specimen after buckling 


cr critical (refers to classical theory) 

edge refers to edge of buckle 

el-b elastic bending 

f final 

i initial after loading (before buckling) 

o initial before loading 

out refers to work done by outside atmosphere 

pi plastic 

r remainder 

y yield 

BUCKLING ANALYSIS 


Assumptions Regarding the Buckling Process 

It is assumed that the thin-walled spherical test specimen is en- 
closed in a larger spherical container as shown in Fig. 1. The space 
between the two spheres is filled with water and at the beginning of 
the test the air outside the container, the water between the spheres 
and the air inside the test specimen are at the same atmospheric pres- 
sure p a • In the first phase of the test the pressure of the water is 
increased to p a + Pi through the use of a pump provided with a no- 
return valve. The increase in the pressure expands the large sphere 
and compresses the small sphere; the effective pressure on the test 
specimen is therefore pi - p c if the notation of Fig. 2 is employed. 
As the pressure in the water should be increased very slowly, it can 
be assumed that the compression inside the test specimen takes place 
isothermally. Thus the pressure inside the specimen before and after 
the compression is related to the volume inside the specimen before 
and after the compression through the equation 

p v = (p + p )v, (l) 

^a o ^a c i 

if the notation of Fig. 2b is used. 
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Under the conditions just described, the specimen is assumed to be 
in a neutral equilibrium. It suddenly jumps over into a buckled shape 
characterized by a dimple of a semivertex angle CC (see Fig. l). The 
dimple is again spherical and it has the same radius as the original 
spherical specimen would have if it were subjected to the prevailing 
pressure difference. The deformations during the buckling process con- 
sist therefore of a uniform expansion corresponding to the change in the 
pressure difference and an inextensional reversal of the curvature of 
the buckled portion of the median surface of the spherical specimen. 

After buckling the water pressure is p a + pf (see Fig. 2c), the 
pressure inside the specimen is p a + p cc , and the effective pressure 
on the spherical specimen is p~ - p • The volume between the two 
spheres does not change during buckling as the water is incompressible. 
The change in the pressure inside the specimen is governed by the 
adiabatic law 


, y 

p a + P cc \\ 

P a + P c ' Vf) 


where for air the value of y is 

y = 1.4 


( 2 ) 


(2a) 


Strain Energy in Unbuckled State 

If the letters R and r denote the median surface radii of the 
outer and inner shells, and the letters H and h the uniform wall 
thicknesses of the outer and inner shells, the uniform normal stresses 
in the shell walls are 

a = p^R/SH in the outer shell 

cr = (p ± - P c )r/2h in the inner shell (3) 

The elastic strain energy stored in the outer shell is 
U = (ct 2 /2E)4ttR 2 H = (2rm 2 H/E)(p i R/2H) 2 = (V 2 )p^VeH 


The total energy stored initially in the two shells is therefore 

I Pi 7T (Pi - Pc) 2 r^ 

1 2 E H 2 E h 


00 
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Elastic Membrane Strain Energy after Buckling 

After buckling the final pressure p a + p f prevails between the 
two shells. Even though part of the inner shell is now concave, the 
absolute magnitude of the membrane stress is the same throughout the 
inner shell as the magnitude of the curvature is the same everywhere. 
The total membrane strain energy stored in the system is 

TT £f TT ( p f ~ p cc } r^ ( 5 ) 

U f ~ 2 E H 2 E h 


Elastic Energy of Bending of Cap 

After buckling the portion of the shell, called the cap, that jump- 
ed over into a position of opposite curvature, is under the action of 
bending moments. The magnitude of the elastic strain energy stored 
because of bending can be evaluated approximately if the energy necessary 
to bend a flat circular plate of a radius r sin “ into a spherical cap 
of radius r is calculated and multiplied by k. The curvature k caused 
by uniform bending moment resultants N distributed along the entire 
circular edge of the plate and acting perpendicularly to the edge is 


r = [12(1 - v)/Eh 3 ]N 


(6) 


Since the curvature sought is l/r, Eq. 6 can be solved for 
r = l/r can be substituted in the result obtained: 


N = 


Eh 3 / r 

12(1 - v) 


N and 


( 7 ) 


Since the edge of the plate must rotate through the angle 2« , and 

since the work done by the edge moment resultant 2N is (l/2)2N2<2 per 
unit length of the circumference, the total work W done is 


W = 2[l2(l - ^)] _1 (Eh 3 / r) a 2ir r sin a 

As the work done is equal to the energy stored, the fol] owing expression 
is obtained for the energy of elastic bending: 


U 


irE 


el .b 3(1 - v) 


h 3 a sin a 


( 8 ) 
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Plastic Work Performed 


In agreement with the assumptions the edge of the buckled cap de- 
forms as a yield hinge. When a yield hinge forms in a ideally rigid- 
plastic material, half the thickness of the sheet is subjected to the 
tensile yield stress, and the other half to the compressive yield stress 
as indicated in Fig. 3. If the two yield stresses are equal numerically, 
the bending moment per unit length is 

M = a (h/2) 2 (9) 

pi y 

This moment is constant during the deformations because no plastic de- 
formations can take place in a rigid-plastic solid until the bending 
moment reaches the value given in Eq.* 7* Hence the total plastic work 
done amounts to 

W - = M „ 2a 2tt r sin a 
pi 

where 2a is the total angle of rotation and 2 ttt sin a is the length 
of the yield hinge. Substitution yields 

2 

W = 7 TO h r a sin a (10) 

pi y 


Changes in Volume of the Two Shells 
Hie volume between the two shells in their natural state is 

V Q = (4tt/3)(R 3 - r 3 ) (n) 

After the water between the two shells is brought to the pressure 
p a + p^ , the outer shell expands and the inner shell is compressed 
with the result that the volume increases. Since the strains in the 
outer and inner shells are 

e = (1 - v)(o/E) = [(l - v)/2](p R/EH) in the outer shell 
e = -(l - v)(cr/E) = -[(l - v)/2](p - p c )(r/Eh) in the inner shell (12) 

the volume of air just before buckling is 
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(13) 


utt r 3 r i-v *£? 3 r i^v (p r p ch 5 ^ 

Vi - 3 ( R [ 1+ 2 ' r L 2 m J i 

me volume V- after buckling can be calculated in a similar manner 
St in addition the volume of the buckle itself must also be taken into 

account . 


Die volume of the cap is 

V c = (l/3)Trr 3 (cos 3 a - 3 cos a + 2) 

Hence the volume V f after buckling is 



1 + 


L-V 


2 EH 


-t . _ 

+ (2rr/3)r (cos a - 3 cos a + 2) 


1 - 


l-v (Pf-P cc )r 


Eh 


me volume of the test sphere immediately before buckling is 

tar 3 I”. l-v «i-V r ! 3 

V, = 


tar 3 r i-v ( v p <> T 

\ = T r I “2 Eh J 


and after buckling 


tar sT, l-v (p f- p C e )r ~| 3 
' 3 r 2 Eh J 

- (27f/3)r 3 (cos 3 a - 3 cos a + 2) 


W 


( 15 ) 


(16) 


( 17 ) 


Work Done by the Air 

me air outside the container remains at the atmospheric pressure 
during the buckling process. Hence the work done by the outside air is 


Vt = t r3 p a{£ + ¥ tf] -[^Tr]} 


(18) 


me work necessary to compress adiabatically the air inside the test 
specimen is 


- 8 - 



( 19 ) 


w 


ad 


<v p c )v i 


y 




+ P cc\ 
+ p c/ 


Ili-1 

y 


Elastic Strain Energy in Edge Zone 

It is easy to see that the deformations of the cap and the remain- 
der in the state of membrane stress existing after the buckle was formed 
are not compatible. The remainder is in compression and the cap is in 
tension. The relative displacement A between cap and remainder can be 
easily calculated with the aid of Figs. 1 and 2: 


A = (1 - V)(r 2 /h) sin a[(p f - P cc )/E] ( 20 ) 

In order to reestablish continuity, edge forces P (perpendicular to 
the axis of the cap) and edge moments M must develop. These can be 
calculated with the aid of formulas derived in Reference 6. The dis- 
placement d perpendicular to the axis of the cap and the rotation 
of the end tangent are given by 


d = 


k, M M 
dM 


+ k dP P 


X “ k XM M + k XP P 


( 21 ) 


with the sign convention shown in Fig. For the cap, the influence 

coefficients can be given as 

k dM = -kyp = 2[3(l-v 2 )]^ (l/E)(r/h 2 ) sin a 

k dP = - 2 [ 3 ( 1 - v2 )^ 1 ^ (l/E)(r/h) 3 ^ 2 sin 2 a 

k *M " M3(1-V 2 )] 3A (l/E)(r l/2 /h 5/2 )C 22 (22) 

where the factors C are defined as 
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(23) 


'11 


11 1 + 2R* 


- (K 


1/2 


'22 v 22 


1 + 2R 


C 12 K ' 


12 1 + 2R* 


The coeficients K were derived by Baker and Cline 7 ; they are tabulated 
in the reference quoted.^ It is of interest to note that they vary 
little with a but change more rapidly with the quantity 

x = q Q! (2^0 

Where 

q =/2 [3(1 - V 2 )] lA (a/b) l/2 (25) 

As the ratio h/a approaches zero, the values of Kqi , K22 > an< i %12 
approach the values n/2 , J~2 , and 1 , respectively. Moreover 

R* = a/a (26) 

cr 


with the classical critical stress 

a = [3(1 - v 2 )]- 1 / 2 E(h/a) (27) 

cr 

The influence coefficients for the remainder are given by the same 
expressions except for the sign; In the expression for k^p the right- 
hand member is positive and in the expression for k^ the right-hand 
member is negative. Equations 23 remain unchanged but the values of the 
Baker-Cline coefficients must be taken from the table established for the 
remainder, and not for the cap. The value of R* also changes sign; it 
is positive for the cap where the stress is tensile and negative for the 
remainder where the stress is compressive. 

With the sign convention used 

P = P = P M = -M = -M (28) 

c r c r 

One condition of compatible deformations is that the lack of fit A 
must be equal in magnitude to the difference in the radial edge dis- 
placement d of the remainder and the cap: 
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( 29 ) 


A = d - d 
r c 


Finally, the sign convention for 


that 


X 

c 


the rotation of the edge tangent is such 

+ X =0 (30) 

r 


Solution of these equations yields 
- 1/2 


p % -1 /p , 2 / ^I2c + ^12r 

M - -(l/2)[3(l-v )] E<h /r)cosec a 2(c 4 C )(c + C llr )-(C 12c+ C 12r )=! 


[3(l-v 2 )] -1 ^ E{h/r)X 2 C osec 2 a 


C 22c + C 22r 


2 ( C 22c +C 22r )(C llc +C llr ) - (C 12c +C X2r ) ‘ 


(31) 


Since the work done by the edge moment resultants is zero as the 
relative rotation between the edges of cap and remainder is zero, work 
is done only by the edge force resultants P • This work is 

W , = (l/2)2irr sin a P A 
el 

This work is equal to the elastic 
stitutions yield 

(l-v) 2 (p f- p cc )2 
T « 

(32) 

It is of interest to note that the value of the last fraction approaches 
1/2 as the final pressure pf - p cc (and thus R* ) decreases and as 
the ratio h/r decreases and approaches zero. 


energy stored in the edge zone. Sub- 


7/2 


sin a 2 

! 2 ( C P? f ' +C ??T'^(^1c + Slr^ C la/’ G 12r'' 


C 22c + C 22r 


Buckling Condition 

Obviously buckling cannot occur unless the energy stored in the 
system before buckling plus the potential of the external loads is at 
least as large as the energy stored in the system after buckling plus 
the work done in an irreversible manner and transformed into heat. Hence 
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a lower limit on the buckling pressure 
energy equation 


U. + 

1 


w 


out 


U f + U el.b 


P i can 


+ U _ + 

edge 


be obtained by solving the 


W + W _ 

pi ad 


(33) 


In addition the geometric constraint 


V 


f 


V. 

1 


(34) 


must be satisfied and the conditions connecting the changes in volume of 
the air contained in the test specimen with the changes in pressure 
(Eqs. 1 and 2) must be observed. The equations contain the nine un- 
knowns p i , Vf > P c > P cc > v i > v f , Vj_ , Vf and a . For their 
solution there are available four equations of geometry (Eqs. 13* 15 > l6, 
and 17); two equations of change of state of the air (Eqs. 1 and 2) , one 
equation of incompressibility of the water (Eq. 3*0 and one energy equa- 
tion (Eq. 33). If a value of a is arbitrarily selected, the eight 
equations suffice to determine all the remaining unknowns. One can plot 
then p^ as a function of a and consider the minimal value of Pj_ as 
the lower bound on the buckling pressure. 


Numerical Values 

The calculations outlined in the preceding sections were programmed 
for the Burroughs 220 digital computer in order to obtain some numerical 
results. For prescribed values of ot y the semi vertex angle of the 
buckle, the gauge pressure p^ in the container was calculated and 
plotted. Figure 5 contains six such curves corresponding to different 
ratios of the radius R of the container to the radius r of the 
specimen. The six values are 1.1 , 1.2 , 1.3 > 1*6 , 2.0 , and 

2.2. The other parameters of the problem are fixed in the following 
manner: 


E = 15 x 10 6 psi a = 30,000 psi 

y 

h/H =0.1 r/h = 500 v = 0.3 


It is to be expected that buckling will occur at the minimal value 
of the gauge pressure. Hence the test arrangement corresponding to 
R/r = 1.2 yields a buckling pressure p^^ which is about 26.5 percent 
of the classical value. This is in good agreement with the test results 
obtained by Sechler and Bollay with the brass shell mentioned in the 
Introduction. The semi vertex angle is about 11 degrees according to 
Fig. 5. 
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As the volume enclosed between test specimen and container is in- 
creased, the semivertex angle increases slowly and the pressure ratio 
decreases slowly. For R/r =2.2 the values obtained are a = 13 degrees 
and Vj_/v cr = 0 . 18 . 

The effect of a change in the material of the specimen (and of the 
container) is illustrated in Fig. 6 . When R/r = 1.2, the minimal 
values of the gauge pressure of buckling amount to approximately 25 , 26 
and kl percent of the classical values of the critical pressure for 
carbon steel, brass and an aluminum alloy, respectively. The correspond- 
ing semivertex angles are about 11 , 11 , and 12.5 degrees. 

To illustrate details of the buckling process, the following data 
are given: 

Characteristics of brass test specimen and container: 

R/r = 1.2 h/H = 0.1 r/h = 500 

E = 15 x 10 ^ psi a = 30,000 psi 


Semivertex angle of buckle 

Gauge pressure before buckling 

Gauge pressure after buckling 

Gauge pressure inside specimen 
after buckling 


a = 11 degrees 
= 19*30 psi 
P f = 7*39 psi 

p cc = 0.012 psi 


The relative importance of the various energy quantities can be seen 
from the following data: 


(l/mr 3 ) (lb - U ) =6.392 x 10 -3 in. lb per in. 3 

(l/7nr J )U , = 2.170 x 10 -3 in. lb per in. 3 
' el.b 

-j /To 

(l/TTYpU^ = 2.36l x 10 " in. lb per in. 

(l/TTr 3 )W ^ = 1 . 69 k x 10 -3 in. lb per in. 3 

o o -5 

(l/mrjW^ = k. 320 x 10 J in. lb per in. 

(l/Trr 3 )W = I .676 x 10 " 3 in. lb per in . 3 
act 
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Simplified Solution 


If use is made of the observation that the changes in volume caused 
by the elastic strains are small compared to the initial volumes , and 
that the semivertex angle OC of the bulge is likely to be small, and 
if in addition the container is assumed to be rigid 


H = oo 

the expressions presented simplify considerably. In particular 

cos 3 a - 3 cos a + 2 = (3/10 
and the volumes of interest become 


Va = 


From Eq. 1 


*a 4tt 3 

v. = — r 

l p + p 3 


(35) 


(36) 


4tt 

3 

r- 

. 3(1 

- v) 

(p. - P 
1 c 

)r- 



(37) 

3 

r 

1 

2 

Eh 

' 111 1 



4ir 

r 3 

r 3(1 

- v) 

(p f - p =c 

)r i 

TT 3 

~4 

(38) 

3 

X 

_L “ 

2 

Eh 

j 

2 r 

u 

4t 

3 

(R 5 

3x 

- r ) + 

2 rr(l 

" v)( Pi - 

) rt 

P c' Eh 


(39) 

4tt 

3 

(R 3 

3\ 

- r ) + 

2tt( 1 

- v)(p f - 

P cc 

4 

is* 

fr 3 ^ 

(40) 


(4i) 


Equation 2 can be re-vritten as 


P + P 
a c 

P + P 
a ^cc 


v. 

l 


and the incompressibility condition becomes 
P i - p c - Pf - p oc + g fl 7 v) T ^ 


(42) 


(43) 


These equations suffice ho express all the quantities of interest in 
terms of the geometric and physical constants of the problem and of the 
initial pressure p i and the angle a . Consideration of the energy 
equation yields then pj_ as a function of a . 
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The energy terms become: 


(p. - P ) 4 

_ 7T 1 c r 

1 " 2 E h 


IT ^ P f ” P cc^ r^ 


U f " 2 E 


U 


7TE 


el.b 3(1 - v) 


,3 2 

h a 


2 2 

W , = iro h ra 

pi y 


w 


ad 


(p + p )v. 
'-^a c i 

7 - 1 


M ' Y - 1 


W = 0 
out 


(44) 

(45) 

(46) 

(47) 

(48) 

(49) 


Moreover, in a first approximation, Uedge can I 16 disregarded. Two 
immediate consequences of the incompressibility condition and Eq. 42 are 


v. 

i 


P„ = P. 


cc 


(50) 


and it follows from Eq. 48 that 

W = 0 (51) 

ad 


From Eqs. 44 and 45 it follows that 

U i - U f = (ir/2)(rVEh)(P i - P f )(P i + P f “ 2P C ) 


But the condition 


V. 

l 


V 


f 


implies 

2jr(r7Eh)(l - v)(p. - p f ) = (Tr/2)r (53) 

while the first of Eqs. 50 leads to the same requirement. Hence, Eq. 33 
becomes simply 
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( 54 ) 


u. 

i 


u =u, . + w 

f el.b pi 


or 


T T r_ / 

2 Eh P i 


P f )(P i + P f 


o v TrE 

p c } 3(l - v) 


h 3 */ 


+ 7 r a h ra (55) 

y 


Eq. 55 together with Eq. 53> can be solved for (p^ - P c ) : 


4(1 - v) 


P i ‘ P c = „2 


E 

3(1 - v) 


(iK(b 


E h 4 
8(1 - v) r 


+ 7TTT — r - a (56) 


From Eqs. 4l and 37 one obtains 


P 0 (p. - p )r 

— = i _ 3 (1 . v) _i 

p + p 2 V ; Eh 


(57) 


a c 


that i s 


1 + (p c /p a } 


j, (1 -v) 


(P, 


P c )r 


Eh 


or approximately, 


P c ? ( p i " p c )r 

_£ = A (1 . v — — 

p 2 V ' Eh 


that is 


P 


c 


- (3/2) (1 - v)(r/Eh) ^ 

(1/pJ + (3/2)(l - v)(r/Eh'5 P i 

a 


Substitution in Eq. 56 yields 


(58) 


p i 


r n 

! ,(3/2) (1 - V)(r/Eh) 

4(1 - v) 

E A\ 3 +o A\ 2 

1 " (l/P ) + (3/2) (1- v) (r/Eh) 
a 

2 

a 

3(1 -V) \rj y\r J 


+ 


E 

8(i rvj 


h J- 

— a 

r 


or 
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u 


l+(3/gU-V)(r/Eh)p, 


1+ ( 1- v ) f E /h\ 

J a 2 (r) 


+ a 


(h) 2 "L E h 4 

AtJ JonTaJ - a 


(59) 


This can also be written as 

r 


Pr Ml-v) 

Pp i - — ~ 

a 


E /h\ 3 

3ii^y w a y 


©•) 


E h 4 




- a 


(60) 


where f3 is the expression in brackets multiplying p^ in Eq. 59; its 
value is almost exactly equal to unity. To obtain Pi m ^ n se t d.Pj_/da 
equal to zero and. solve for a : 


6 16(1 - v) g h 

m E r 


E 


-7 v- + a 

3(l - v) r y 


(61) 


Substitution in Eq. 60 yields 


Pi 


min 


I 7T~ — r (Vp) 

o (1 - v)r mm ' 


(62) 


A numerical example will show the accuracy of the approximate formula. 
If a brass sphere is characterized by the values 

,6 


E = 15 x 10 


a = 30,000 

y 


h/r = l/ 500 

v = 0.3 


one obtains from Eq. 6l 


a 6 = 46.2 X 10" 6 a = 0.189 = 10° 45’ 


and from Eq. 62 


p = 20.2 psi P i /P cr = (20.2/72.3)100 = 27.9$ 


SIMPLE BUCKLING FORMULAS 


If the assumptions of the last section are maintained, namely that 
the container is perfectly rigid, and in addition v is taken as 0.3; the 
results of the analysis can be presented in the following simple formulas: 


a . 

min 


1^7 


cr 


^ 6 /a 

[ cr 


/ v 


+ 1.27 



(63) 
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(64) 


a . = l.2k f 

min V 

p • . = k .16 a 


hV /6 A + 2 . 

r 


1/6 


■ imin 


r ) e ) 

°„ r \2/3 /«„. «„\ 2 /3 


cr \ 
cr E J 


cr 

E 


4 - 1.27 


0 \ 


= 1-26 E( + 2.1 /3 


■ imin 


1 min 


= 0.630 + 2.1 ^ ^ 


(65) 

( 66 ) 

(67) 


where cr is the critical stress defined in Eq. 27* 
cr ^ 

The last equation shows that for comparatively thick shells with a 
low value of the yield stress the buckling stress decreases more rapidly 
with decreasing h/r ratio than is predicted by the classical small- 
deflection formula, while for comparatively thin shells with a high 
value of the yield stress the opposite is true. 
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FIG. I. BUCKLED SPHERICAL TEST SPECIMEN 

IN SPHERICAL CONTAINER 
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FIG. 2 THREE PHASES OF BUCKLING 



FIG. 3. STRESS DISTRIBUTION IN WALL OF 
SHELL CORRESPONDING TO PLASTIC 
HINGE. 


- 22 


(a) Buckled cap 





(b) Remainder 


« 



FIG. 4. CONVENTION FOR EDGE DISPLACEMENTS 
AND STRESS RESULTANTS. 
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FIG. 5 EFFECT OF VOLUME OF CONTAINER ON BUCKLING PRESSURE 





FIG. 6. EFFECT OF MATERIAL OF SPECIMEN ON BUCKLING PRESSURE 




